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HIGHLIGHTS 


►  We  develop  a  ID  mathematical  model  of  lithium-ion  battery. 

►  We  use  electrical  analogy  to  solve  mass  and  charge  transport  equations. 

►  The  model  is  implemented  in  a  standard  software  used  in  electrical  engineering. 

►  The  model  can  be  easily  employed  for  the  simulation  of  electrical  power  systems. 
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This  article  deals  with  an  electrical  model  of  lithium-ion  battery,  based  on  a  ID  analog  representation  of 
mass  and  charge  transport  phenomena.  By  using  an  electrical  analogy  to  describe  migration  and  diffu¬ 
sion  of  species  and  charges  in  electrodes  and  electrolyte,  the  model  can  be  directly  implemented  in 
standard  simulation  software  used  in  electrical  engineering  (such  as  Saber,  in  our  case),  so  that  it  can  be 
easily  employed  for  the  simulation  of  electrical  systems  involving  lithium-ion  batteries.  The  paper 
explains  how  the  analog  model  is  obtained  from  mass  and  charge  transport  equations,  then  it  presents 
a  comparison  with  results  obtained  by  means  of  a  finite  element  method.  At  last,  some  electrical  power 
systems  are  simulated  in  order  to  show  advantages  of  our  modeling  method. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium-ion  batteries  are  based  on  insertion  electrodes,  so  that 
their  electrical  behavior  is  marked  by  line  effects  due  to  insertion  or 
deinsertion  of  lithium  in  electrode  active  material.  As  demon¬ 
strated  in  Ref.  [1],  the  slow  kinetics  of  these  phenomena  limits 
performances  of  lithium-ion  batteries,  particularly  in  terms  of  quick 
charge.  It  can  also  lead  to  unexpected  effects,  such  as  plating  (metal 
lithium  deposition  at  electrode  surfaces).  It  is  therefore  crucial  to 
develop  a  lithium-ion  model  that  combines  a  local  description  of 
physical  phenomena  with  an  implementation  suitable  for  simula¬ 
tion  of  electrical  power  systems. 

Within  a  lithium-ion  cell  take  place  key  processes  that  are 
strongly  coupled:  lithium  diffusion  in  electrode  active  material, 
lithium  ion  migration  and  diffusion  in  electrolyte,  electron 
conduction  in  electrodes,  and  electrochemical  kinetics  at  electrode/ 
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electrolyte  interface.  Many  mathematical  models  have  already  been 
developed  in  order  to  describe  locally  these  phenomena,  through 
coupled  partial  differential  equations  (PDE)  involving  space  and 
time.  One  of  the  earlier  studies,  by  Doyle  et  al.  [2],  led  to  an 
isothermal  and  one-dimensional  model,  relative  to  a  metal  lithium 
cell.  In  another  significant  early  work  [3],  the  same  authors 
extended  the  approach  to  the  modeling  of  dual  lithium  ion  inser¬ 
tion  cell.  These  works,  on  which  our  1 D  lithium-ion  battery  model 
is  based,  have  been  recently  improved  by  Smith  and  Wang  [4]  and 
by  Cai  and  White  [5],  who  introduced  thermal  heat  transfer, 
temperature  dependant  parameters  and  electrothermal  coupling. 
According  to  results  presented  in  the  literature,  mathematical 
lithium-ion  battery  models  are  obviously  accurate,  and  enable  to 
study  specific  local  effects.  However,  they  are  hardly  usable  in 
a  system  approach,  for  which  a  lithium-ion  battery  model  easy  to 
handle  is  required. 

This  is  the  reason  why  many  electrical  circuit  models  have  been 
developed,  most  often  based  on  the  Randles  equivalent  scheme 
that  represents  the  electrical  behavior  of  an  electrochemical 
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electrode-electrolyte  interface.  This  scheme  is  generally  composed 
of  a  series  resistance  for  purely  ohmic  voltage  drops  in  the  elec¬ 
trode  and  the  electrolyte,  an  electric  double  layer  capacitance  in 
parallel  with  faradic  impedances  that  describe  charges  transfer 
phenomena  that  occur  at  the  electrode-electrolyte  interface,  and 
active  material  transport  in  electrode  and  electrolyte  [6,7].  In 
electrical  engineering,  this  basic  Randles  equivalent  scheme  is  often 
reduced  to  a  constant  voltage  source  in  series  with  a  constant 
resistance.  Sometimes,  like  in  Ref.  [8],  a  parallel  RC  dipole  is  added 
in  series,  to  describe  the  electric  double  layer  phenomenon. 
However,  many  models  improve  this  basic  representation  by 
implementing  parameters  depending  on  state  of  charge  (SOC).  This 
is  for  example  the  case  in  Ref.  [9],  in  which  a  linear  variation  versus 
SOC  is  considered  for  both  open-circuit  voltage  (OCV)  and  series 
resistance.  In  Ref.  [10],  battery  voltage  is  computed  by  a  third-order 
polynomials  versus  SOC,  for  different  discharge  current  values.  An 
exponential  empirical  law  is  used  by  Tremblay  et  al.  in  Ref.  [11  ],  to 
describe  OCV  variation  versus  SOC.  The  resulting  circuit  model  has 
been  implemented  in  Refs.  [12]  and  [13],  to  study  power  generation 
systems  using  lithium-ion  storage.  This  approach  has  been 
improved  by  Tofighi  and  Kalantar  [14],  who  extended  the  expo¬ 
nential  dependance  versus  SOC  to  series  resistance  and  short  and 
long  transient  parallel  RC  dipoles.  In  Ref.  [15]  at  last,  Urbain  et  al. 
developed  an  accurate  energetical  modeling  of  lithium-ion 
batteries  including  SOC  and  current  dependances,  electrode 
porosity  and  concentration  effects,  and  relaxation.  This  kind  of 
models  can  give  accurate  results,  but  they  are  most  often  empirical, 
and  of  course  can  not  describe  local  effects. 

Another  way  to  handle  the  Randles  equivalent  scheme  is  to 
develop  dynamic  models  based  on  harmonic  response.  This 
approach  is  used  in  Refs.  [16,17]  for  modeling  lead-acid  batteries,  in 
Ref.  [18]  for  modeling  NiMH  batteries,  in  Ref.  [19]  for  modeling 
lithium-ion  batteries,  and  in  Ref.  [20]  for  modeling  the  electrical 
behavior  of  a  PEM  fuel  cell.  Thus,  this  approach  is  widespread  in 
electrochemical  device  modeling.  Nevertheless,  although  it  is  most 
often  physics-based,  it  fundamentally  generates  only  small  signal 
models,  which  are  theoretically  not  available  in  large  signal 
conditions,  because  of  non-linearities  of  electrochemical  devices. 
Moreover,  in  the  case  of  batteries,  it  is  difficult  to  obtain  an 
equivalent  circuit  at  a  mean  current  different  from  zero,  because  of 
the  SOC  dependence  of  battery  behavior.  And  once  again,  this  kind 
of  models  can  hardly  describe  local  effects.  However,  we  can  find  in 
the  literature  battery  models  in  Laplace  domain  that  are  deduced 
from  a  mathematical  PDE  formulation,  such  as  in  Ref.  [21  ].  In  this 
case,  local  effect  description  can  be  associated  with  system  oriented 
models.  This  is  the  purpose  of  our  work,  but  directly  in  time 
domain. 

Indeed,  we  focus  in  this  article  on  a  lithium-ion  battery  cell 
electrical  model  implemented  in  a  simulation  software  commonly 
used  in  electrical  engineering  to  design  systems.  This  model  is 
isothermal  and  includes  transient  description,  double-layer 
phenomenon,  and  1 D  mathematical  representation  of  charge  and 
mass  transport  phenomena  in  electrodes  and  electrolyte.  For  this 
purpose,  we  use  an  electrical  analogy  of  species  and  charge 
migration  and  diffusion  equations,  so  that  the  model  can  be  directly 
implemented  in  a  simulation  software  dedicated  to  electrical 
circuit.  It  can  be  therefore  easily  employed  for  the  simulation  of 
electrical  systems  involving  lithium-ion  batteries. 

In  the  first  part,  mass  and  charge  transport  equations,  electro¬ 
chemical  kinetic  equations,  and  boundary  conditions  of  our 
lithium-ion  battery  model  are  recalled.  Then,  electrical  analogy  of 
mass  transport  equations  in  electrodes  and  electrolyte,  and  of 
associated  boundary  conditions,  is  detailed,  as  well  as  the  model 
assembly.  A  numerical  validation  is  carried  out,  by  comparing  our 
model  results  to  those  obtained  by  means  of  a  finite  element 


resolution  method.  Finally,  some  applications  are  simulated  (series 
operation,  association  with  a  power  electronics  converter)  in  order 
to  show  advantages  of  our  modeling  method. 

2.  Mathematical  ID  model  of  lithium-ion  batteries 

A  lithium-ion  cell  is  composed  of  two  electrodes,  an  electrolyte, 
a  porous  separator,  and  two  current  collectors.  We  briefly  describe 
hereafter  materials  commonly  used  for  lithium-ion  technology,  and 
principles  of  operation.  Further  information  can  be  found  in  Ref.  [22]. 

Active  material  of  the  positive  electrode  is  generally  a  transition 
metal  oxide  containing  lithium,  mainly  cobalt  based  (LiCoC^).  The 
insertion  of  lithium  operates  under  a  potential  ranging  between 
3.5  V  and  4  V  (compared  to  lithium).  NMC  materials  (LiNiyMnz_ 
Coi_y_z02)  also  penetrate  the  market  of  portable  electronics, 
because  they  enable  to  reduce  the  content  in  expensive  cobalt.  One 
also  uses,  in  particular  for  industrial  applications,  metallic  oxides 
containing  nickel  (LiNi02)  or  iron  phosphate  (LiFeP04,  also  called 
LFP).  The  insertion  material  of  the  negative  electrode  is  a  carbon 
material  of  coke  type,  graphite,  or  amorphous  carbon.  The  insertion 
of  lithium  operates  under  a  potential  ranging  between  0  V  and  0.3  V 
(compared  to  lithium). 

The  electrolyte  ensures  ion  migration  within  the  battery.  In  the 
case  of  lithium-ion  accumulators,  it  is  a  mixture  of  cyclic  and/or 
linear  carbon  solvents,  such  as  propylene  carbonates  (PC),  in  which 
a  lithium  salt,  lithium  hexafluorophosphate  LiPF6  for  example,  is 
dissolved.  Moreover,  the  reduction  of  the  electrolyte  solvents 
enables  to  generate  a  passivation  layer  (known  as  SEI,  for  Solid 
Electrolyte  Interphase)  at  negative  electrode/electrolyte  interface, 
thus  avoiding  the  degradation  of  the  electrolyte  under  too  low 
potentials,  and  protecting  carbon  material. 

The  separator  constitutes  a  physical  barrier  between  the  two 
electrodes.  Thus,  it  must  be  an  electronic  insulator.  It  must  also  be 
porous,  in  order  to  allow  the  ion  migration  between  electrodes. 
Separators  used  in  lithium-ion  technology  usually  contain  poly¬ 
ethylene  (PE)  and/or  polypropylene  (PP). 

Lithium-ion  accumulators  operate  according  to  the  “rocking 
chair”  principle  [22]:  during  operation,  charge  migration  within  the 
device  is  ensured  by  lithium  ions  Li+,  which  leave  the  first  electrode 
structure  to  insert  into  the  second  electrode  structure.  This 
phenomenon  of  insertion-deinsertion  of  lithium  ions  in  electrode 
materials  can  be  described  at  the  negative  electrode  by: 

discharge 

Liy+dvc6  <  -  Liy C6  +  dy  •  Li  +  dy  •  e- ,  (1) 

charge 

and  at  the  positive  electrode  by: 

,  discharge 

LizM02  +  dz-Li+  +  dze-  ^  ^  Liz+dzM02,  (2) 

charge 

y  (0  <  y  <  1)  and  z  (0  <  z  <  1)  being  lithium  insertion  rates  in 
negative  and  positive  electrodes,  respectively. 

In  this  section,  we  first  recall  assumptions,  transport  equations 
and  boundary  conditions  of  a  mathematical  lithium-ion  battery 
model  first  introduced  by  M.  Doyle,  T.  F.  Fuller  and  J.  Newman  in 
1993  for  a  metal  lithium  cell  [2],  in  1994  for  a  dual  lithium  ion 
insertion  cell  [3],  and  recently  improved  by  Smith  and  Wang  [4]. 
One  will  find  in  Ref.  [1]  an  experimental  validation  of  this  mathe¬ 
matical  lithium-ion  battery  model,  carried  out  by  K.  Smith  and  C.-Y. 
Wang  on  a  6  Ah  element.  As  indicated  in  Fig.  1,  the  ID  lithium-ion 
cell  model  is  composed  of  three  domains: 

1)  the  composite  negative  electrode,  which  contains  active 
material  particles  (LiyC6)  impregnated  with  electrolyte, 
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Fig.  1.  1 D  lithium-ion  cell  model. 


Rs  being  the  particle  radius,  assumed  constant  in  each  electrode. 
At  last,  we  can  associate  to  GLl  a  volume  rate  of  Li+  current  gener¬ 
ation  (in  A  m-3)  as  follows: 

jLi(x)  =  as  jo  ■  (exp  («of))  -  exp (  -  ) .  (5) 

In  discharge  regime,  and  contrary  to  charge  regime,/1  is  positive 
at  the  negative  electrode,  and  negative  at  the  positive  electrode. 


2.2.  Solid  phase  equations 


2)  the  separator,  a  porous  dielectric  that  ensures  electronic  insu¬ 
lation  and  ion  migration, 

3)  the  composite  positive  electrode,  which  contains  active  mate¬ 
rial  particles  (LizM02)  impregnated  with  electrolyte. 


2.2.1.  Lithium  transport 

Lithium  transport  in  active  material  particles  constituting  the 
solid  phase  of  the  electrodes  is  assumed  to  be  purely  diffusive.  In 
a  ID  approach,  the  flux  density  of  lithium  in  a  particle  can  be 
written  as: 


Assumptions  and  simplifications  adopted  in  the  present  model 
are  as  follows: 


cps(r)  =  —Ds  ■ 


9cs(r) 

dr 


(6) 


1)  the  cell  temperature  remains  constant  and  homogeneous 
throughout  the  cell, 

2)  in  the  electrolyte  phase,  mass  and  charge  transports  are  one¬ 
dimensional,  along  the  x  direction  in  Fig.  1.  They  are  both  due 
to  diffusion  (Fields  diffusion  law)  and  migration  (Ohm’s  law). 
The  convection  is  neglected, 

3)  in  the  solid  phase,  active  material  particles,  in  or  out  which 
lithium  diffuses,  are  supposed  to  be  spherical,  so  that  mass 
transport  due  to  diffusion  (Fields  diffusion  law)  is  treated 
spherical  in  this  phase.  Particle  radius  is  considered  small 
enough  compared  with  electrode  thicknesses,  so  that  lithium 
concentration  at  solid/electrolyte  interface  is  computed  as 
a  continuous  function  of  the  space  variable  x.  Electron  trans¬ 
port  through  migration  (Ohm’s  law)  is  one-dimensional,  along 
the  x  direction  in  Fig.  1. 


A  balance  of  lithium  in  a  particle  allows  to  determine  the 
equation  describing  the  lithium  transport  in  solid  phase: 


9cs(r) 


_  1  Aid 

6t  “  r2  9r  Us 


9cs(r) 

dr 


(7) 


Boundary  conditions  are  of  Neumann  type.  At  the  particle 
center,  lithium  flux  density  is  equal  to  zero: 


<Ps(  0)  =  0, 


(8) 


and  at  the  solid/electrolyte  interface,  lithium  flux  density  is  linked 
to  electrochemical  reaction  rate: 


<Ps(Rs) 


jLi(x) 

as-F' 


(9) 


The  problem  has  four  unknowns:  the  ionic  concentration  ce  of 
the  electrolyte,  the  lithium  concentration  cs  in  active  material 
particles,  the  electrolyte  phase  potential  0e.  and  the  solid  phase 
potential  </>s.  Electrolyte  phase  unknowns  ce  and  0e  are  available  in 
the  three  domains,  whereas  solid  phase  cs  and  </>s  are  only  available 
in  both  electrode  domains. 


2.1.  Electrochemical  reaction  rate 

The  volume  rate  GLl  of  lithium  ion  generation  at  the  solid/elec¬ 
trolyte  interface,  i.e.  the  number  of  Li+  moles  released  in  electrolyte 
phase  per  time  and  electrode  volume  units,  is  determined  by 
Butler- Volmer  kinetics  law: 

GLi(x)  =  jh  ■  (exp  («0f)  -  exp  (  -  ar^/) ),  (3) 

where  jo  is  the  exchange  current  density,  as  is  the  electrode  active 
surface  per  electrode  volume  unit,  p  is  the  local  interface  over¬ 
voltage,  Uj  is  the  thermal  voltage  (ut  =  FT/F,  R  being  the  gas  constant, 
T  the  temperature,  and  F  the  Faraday’s  constant),  a0  and  aT  are  the 
oxidation  and  reduction  transfer  coefficients,  respectively.  Active 
material  particles  being  supposed  spherical,  and  taking  into  account 
the  active  material  volume  fraction  rs,  the  ratio  as  is  given  by: 


2.2.2.  Charge  transport 

In  solid  phase,  charge  transport  corresponds  to  electronic 
conduction.  Denoting  js  the  associated  electronic  current  density,  </>s 
the  potential  of  the  solid  phase  and  crs  eff  the  effective  electronic 
conductivity  of  the  considered  electrode,  electron  transport  in 
electrodes  is  described  by  Ohm’s  law  and  can  be  written  in  ID, 
according  to  the  “x”  axis  shown  in  Fig.  1 : 


,  M  .  00s(x) 

MX)  —  as,eff'  5 


(10) 


the  effective  electronic  conductivity  being  given,  versus  the  mate¬ 
rial  electronic  conductivity  as  and  the  active  material  volume 
fraction  rs,  by: 

fl^eff  =  (11) 

A  balance  of  electronic  charges  in  electrode  enables  to  express  the 
charge  conservation  equation  in  solid  phase  as  follows,  considering 
that  the  generation  rate  of  electrons  in  the  solid  phase  is  obviously 
equal  to  the  generation  rate  of  Li+  cations  in  the  electrolyte  phase: 


f  -  -j“« 


(12) 


Boundary  conditions  at  electrode/separator  interfaces  are 
characterized  by  a  current  density  equal  to  zero,  because  of  the 
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insulating  nature  of  the  separator.  What  is  written,  according  to  the 
notations  defined  in  Fig.  1 : 


fis(^neg)  —  0 

{JsfLceW  ~  ^pos)  =  0 


(13) 


At  electrode/collector  interfaces,  electronic  current  density  is 
imposed  by  the  external  circuit,  which  can  be  written,  if  Jceii  is  the 
cell  current  and  Aceii  the  electrode  surface: 


|7s(0)  =  Ice  n/Aceii 

\js(L)  =  fCell/A:ell 


(14) 


For  the  electrolyte  ionic  conductivity  cre,  which  depends  on  the 
electrolyte  concentration,  we  will  make  use  of  the  following 
empirical  law  [4]: 

ffe(x)  =  0.001 58  •  ce(x)- exp  (-  5.363  x  10-5  c^4(x)).  (21) 

If  the  electrolyte  potential  0e  is  refered  by  a  lithium  reference 
electrode  at  the  average  electrolyte  concentration  (corresponding 
to  the  OCV  measurement  conditions),  the  effective  diffusional  ionic 
conductivity  <7De,eff  °f  th e  electrolyte  is  given  by: 

^De,eff  —  2  •  Uj  •  (tp  —  1)  •  (7e  eff .  (22) 


2.3.  Electrolyte  phase  equations 
2.3 A.  Ion  transport 

Ion  transport  in  electrolyte  phase  is  a  combination  of  two 
physical  phenomena.  The  first  one  is  due  to  the  migration  caused  by 
local  electric  field,  and  the  other  one  is  ion  diffusion  that  results 
from  concentration  gradients.  One  demonstrates  that  species 
conservation  in  electrolyte  phase  leads  to  the  following  partial 
differential  equation,  involving  the  ionic  concentration  ce,  the 
electrode  porosity  te,  the  Li+  transference  number  tp  and  the 
effective  salt  diffusion  coefficient  in  the  electrolyte  phase  De,eff: 


Lastly,  the  balance  of  ionic  charges  enables  to  express  the  charge 
conservation  equation  in  electrolyte  phase  as  follows: 


(23) 


It  can  be  noticed  that  the  cell  current  density  JCeii,  equal  to  the 
sum  of  electron  current  density  Js  in  solid  phase  and  ion  current 
density Je  in  electrolyte  phase,  has  a  divergence  equal  to  0  (Eqs.  (12) 
and  (23)),  and  therefore  is  constant  throughout  the  cell. 

Electrode/collector  interfaces  are  characterized  by  a  purely 
electronic  current  density,  so  that  ion  current  density  is  equal  to 
zero  at  these  bondaries.  Considering  Eq.  (18),  this  can  be  simplified 
in: 


The  effective  salt  diffusion  coefficient  is  given  by  the  Bruggeman 
relation  [23]: 

Deeff  =  E?-De,  (16) 

where  De,  i.e.  the  salt  diffusion  coefficient  in  the  electrolyte  phase, 
is  related  to  anion  and  cation  diffusion  coefficients: 


(0)  =  0 
(^cell)  ~  0 


(24) 


At  electrode/separator  interfaces,  the  continuity  of  ion  current 
density  is  applied. 


2  •  Dpp-  •  PLj+ 
DPF6-  +  DLi+  ' 


(17) 


Boundary  conditions  at  electrode/collector  interfaces  are  char¬ 
acterized  by  a  mass  flux  equal  to  zero,  what  can  be  written  as: 


1>~° 


(18) 


At  electrode/separator  interfaces,  the  continuity  of  ionic 
concentration  is  applied. 


2.3.2.  Charge  transport 

The  total  ion  current  density  Je  associated  with  ion  transport  in 
the  electrolyte  phase  is  the  sum  of  two  conduction  components 
(anion  and  cation  transport  by  migration  induced  by  potential 
gradient)  and  two  diffusion  components  (anion  and  cation  trans¬ 
port  by  diffusion  due  to  concentration  gradient).  It  can  be  expressed 
versus  the  electrolyte  concentration  ce  and  the  electrolyte  potential 
0e  as  follows: 

JeM  —  ^e,effM’  ^De,eff(^)'  •  (19) 

The  effective  ionic  conductivity  cre,eff  of  the  electrolyte  depends 
on  the  electrolyte  ionic  conductivity  and  the  electrode  porosity, 
according  to  the  Bruggeman  relation  [23]: 

°e,eff  =  (20) 


2.4.  OCV,  electrode  overvoltages  and  cell  voltage 

Electrode  equilibrium  potentials  are  local  variables  that  depend 
on  lithium  insertion  rate  at  the  surface  of  solid  phase  particles,  i.e. 
on  the  ratio  cSte(x)/csmax  between  lithium  concentration  at  r  =  Rs 
and  maximum  lithium  concentration  in  solid  phase.  Curves  plotted 
hereafter  in  Fig.  2  have  been  determined  by  means  of  empirical 
formula  given  in  Ref.  [4].  For  implementing  the  model  in  Comsol 
Multiphysics  or  in  Saber,  we  will  make  use  of  tables  associated  with 
these  curves. 

The  local  electrode  overvoltage  involved  in  the  expression  of  the 
electrochemical  reaction  rate/1  (Eq.  (5))  is  equal  to  the  difference 
between  the  interfacial  voltage  and  interfacial  equilibrium  voltage. 
It  is  therefore  expressed  as  a  function  of  the  potential  of  the  solid 
phase  0s,  the  potential  of  the  electrolyte  phase  0e,  and  the  electrode 
equilibrium  potential  U  as  follows: 

V(x)  =  (0s(x)-0e(x))-L/(x).  (25) 

Once  again,  it  should  be  noticed  here  that  in  this  formulation 
commonly  used  in  the  literature,  the  electrolyte  potential  0e  is 
referred  by  a  lithium  reference  electrode  at  the  average  electrolyte 
concentration,  corresponding  to  the  OCV  measurement  conditions. 
As  a  result,  the  electrode  equilibrium  voltage  U  involved  in  Eq.  (25) 
only  depends  on  the  local  insertion  rate  cSie(x)/csmax,  the  theoretical 
effect  of  the  local  electrolyte  concentration  on  electrochemical 
kinetics  being  included  in  the  expression  the  effective  diffusional 
ionic  conductivity  (JDe.eff- 

Finally,  the  cell  voltage,  which  is  the  potential  difference 
between  the  two  current  collectors,  is  therefore  equal  to  the 
difference  of  solid  phase  potentials  at  x  =  Lceii  and  x  =  0: 
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Kell  =  0sC^cen)  -  0s(O).  (26) 


particle  element.  For  more  accuracy  and  for  convenience  regarding 
boundary  conditions,  a  two-order  circuit  is  implemented,  with  the 
following  component  expressions: 


Ivm  —  - t - 

4it(r  +  Ar/2)2 -Ds 
R  Ar/2 

Ll2  47r(r  + Ar)2-Ds  .  (29) 

Cm  =  5*-((r  +  A  r/2)3-r3) 

CLi2  =  |7T-((r  +  Ar)3-(r  +  Ar/2)3) 


Lithium  diffusion  in  solid  phase  is  then  described  by  associating 
in  series  such  equivalent  circuits.  In  the  analogical  model  presented 
here,  a  space  discretization  in  10  elements  is  used  for  each  elec¬ 
trode  (A r  =  Ks/10).  Fig.  4  depicts  input,  output  and  boundary 
conditions  of  the  resulting  lithium  diffusion  model.  Model  input,  i.e. 
electrochemical  reaction  rate  jL\  is  required  for  setting  boundary 
molar  flow  conditions,  according  to  Eq.  (9)  that  becomes  here: 

Ns(Rs,x)  =  4ttR24^.  (30) 

asr 

Model  output,  the  lithium  insertion  rate  at  the  surface  of  active 
material  particles,  enables  the  calculation  of  the  local  electrode 
equilibrium  potential.  Boundary  conditions,  i.e.  Eqs.  (8)  and  (30)  for 
lithium  molar  flows  at  r  =  0  and  r  =  Rs,  are  described  by  means  of  an 
open  circuit  and  a  current  source,  respectively. 


3.  Lithium-ion  battery  modeling  using  electrical  analogy 


3.1.  Modeling  lithium  transport  in  solid  phase 


Space  discretization  of  Eqs.  (6)  and  (7)  leads  to  the  following 
system,  for  the  description  of  lithium  diffusion  in  active  material 
particles  constituting  the  solid  phase: 


cs(r  +  Ar)  =  cs(r)-^^--Ns(r) 

Ns(r  +  Ar)  =  Ns(r)  -  Qjr((r  +  Ar)3-r3)) 


(27) 


Ar  being  the  space  discretization  step,  and  Ns  being  the  lithium 
molar  flow  within  active  material  particles,  defined  in  spherical 
coordinates  versus  the  lithium  molar  flow  density  (ps  as  follows: 


3.2.  Modeling  charge  transports  and  mass  transport  in  electrolyte 
phase 


Electronic  and  ionic  currents  being  denoted  Js  (/s  =  Ace\\  Js)  and  Je 
(/e  =  AreiiJe).  respectively,  space  discretization  of  Eqs.  ((10),  (12), 
(15),  (19)  and  (23)  leads  to: 


/</>s(x  +  Ax)  =  0s(x)-Ks,eirfs(x) 

U(x  +  Ax).  =  fs(x)-/e  ltW 

j  ce(x  + Ax)  =  ce(x)  -  Rjon  ■  Ne(x) 

{  Ne(X  +  Ax)  =  Ne(x)  +  Nion(x)  -  Cion 

r  Mx  +  Ax)  =  4>e(x)  -  Ue,eif/e(x)  -  VeltW 
\  Ie(X  +  Ax)  =  /e(X)  +  4lt(x) 


(31) 


Ns(r)  =  4tt  r2-cps(r).  (28) 

This  non-linear  system  can  be  computed  using  an  electrical 
analogy,  for  which  lithium  concentration  and  molar  flow  are 
assimilated  to  voltage  and  current,  respectively.  Fig.  3  presents  the 
equivalent  electrical  scheme  of  lithium  diffusion  in  a  spherical 


c,(r) 


A 

< - 

RLil  RLi2 

< - 

A 

^Lil  ^Li2 

cs(r+Ar) 


jLi(x) 


Fig.  3.  Equivalent  circuit  of  lithium  diffusion  in  an  element  of  solid  particle. 


Fig.  4.  Electrical  analogy  of  lithium  diffusion  in  electrode  solid  phase  particles. 
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with  the  following  expressions  for  charge  transport  parameters: 


^s,elt  — 


Ax 


Azell  ’  ^s,eff 
R  =  — 

e’elt  A:eirffe,eff  • 

4lt(*)  =  ^cel  l-Ax-jLi(x) 


(32) 


and  the  following  parameters  for  mass  transport  in  the  electrolyte 
phase: 


^ion  —  a 
A 


Ax 


LceH'^e,eff 


Qon  —  £e  '  Azell '  ^ 

^ion(^)  =  Azell' 0  —  ^p) 


(33) 


jLi(x) 


equations  of  Eq.  (31)  and  associated  equivalent  scheme  are  of 
course  not  to  be  considered. 

The  lithium  insertion  rate  cSie(x)/cs max  at  the  surface  of  active 
material  particles,  determined  by  the  model  of  lithium  diffusion  in 
solid  phase,  enters  the  cartesian  model  for  calculating  the  local 
electrode  equilibrium  potential  (by  means  of  look-up  tables  asso¬ 
ciated  with  curves  of  Fig.  2).  Then,  the  local  electrode  overvoltage  p 
and  the  electrochemical  reaction  rate  jLl  are  evaluated  according  to 
Eq.  (25)  and  to  Eq.  (5),  respectively.  This  last  variable  is  required  to 
compute  current  creation  terms  (Je it  of  Eq.  (32)  and  Njon  of  Eq.  (33)), 
and  is  also  used  as  an  input  of  the  solid  phase  diffusion  model. 

Charge  transport  in  solid  phase,  and  mass  and  charge  transports 
in  electrolyte  phase,  are  described  by  associating  in  series  such 
equivalent  circuits.  In  our  analogical  model,  a  space  discretization 
in  10  elements  is  used  for  each  electrode  (Ax  =  <5neg/10  or  (Ax  =  <5pos/ 
10),  and  for  separator  (Ax  =  <5sep/10). 


3.3.  Lithium-ion  battery  model  implementation 


Ax  being  the  space  discretization  step.  This  equation  system, 
available  in  both  electrodes,  can  be  computed  using  the  equivalent 
electrical  scheme  depicted  in  Fig.  5.  In  separator,  solid  phase 


c«(x>/csm,lx 


jLl(x) 


Fig.  5.  Equivalent  circuit  of  charge  and  mass  transports  in  an  element  of  cell  (electrode 
or  separator). 


Fig.  6  presents  the  complete  lithium-ion  battery  model  that  has 
been  implemented  in  Saber  software.  It  combines: 

1)  analogical  transport  models  previously  described,  which 
calculate  concentrations  ce  and  cs,  electrical  potentials  </>e  and 
</>s,  mass  flows  Ne  and  Ns,  electrical  currents  Je  and  Js,  over¬ 
voltage  7]  and  electrochemical  reaction  rate  jLl, 

2)  a  controlled  voltage  source,  which  couples  the  electrical  model 
to  the  transport  model  via  the  voltage  across  the  cell. 

As  depicted  in  Fig.  6,  boundary  conditions  on  electronic  current 
Is  at  electrode/separator  interfaces  (Eq.  (13)),  and  on  ce  and  /e  at 
electrode/collector  interfaces  current  (Eqs.  (18)  and  (24))  are 
treated  by  means  of  floating  points  (corresponding  circuits  are 
opened  at  these  interfaces).  Boundary  conditions  on  Js  at  electrode/ 
collector  interfaces  (Eq.  (14))  are  implemented  by  current  sources, 
equal  to  /cen.  At  last,  continuity  conditions  for  ce  and  Je  at  electrode/ 
separator  interfaces  are  naturally  modeled  by  connecting  corre¬ 
sponding  electrode  and  separator  circuits. 

To  compute  the  transport  model,  a  potential  reference  and 
initial  conditions  on  concentrations  are  needed.  The  potential 
reference  is  set  at  the  negative  terminal  (i.e.  </>s(0)  =  0).  Initial 
electrolyte  concentration  is  chosen  constant  throughout  the  whole 
electrolyte  domain,  and  equal  to  the  average  electrolyte  concen¬ 
tration  ce,mean-  Initial  lithium  concentrations  in  solid  phases  are  also 
chosen  constant  throughout  considered  domains,  and  equal  to: 

cs,init  =  Cs  o  H  ^qq  (cs,100  —  cs,o)5  (34) 

where  SOQnit  is  the  initial  state  of  charge  (in  %),  and  cSio  and  cSiioo 
are  the  lithium  concentrations  in  solid  phase  when  the  cell  is  fully 
discharged  or  charged,  respectively. 

4.  Comparison  with  the  model  solved  by  finite  element 
method 

The  ID  mathematical  model  of  lithium-ion  battery  described  in 
section  2  has  also  been  implemented  in  commercial  software  that 
enables  multiphysics  modeling  (Comsol  Multiphysics,  in  our  case). 
The  resolution  is  based  on  the  finite  element  method,  each  domain 
is  discretized  in  40  elements.  The  aim  of  this  section  is  to  compare 
results  given  by  this  finite  element  multiphysics  software,  to  those 
obtained  by  means  of  a  analogical  representation  of  transport 
phenamena  implemented  in  Saber  software. 
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Icell(t) 


Table  1 

Model  parameters  [4]. 


Parameter 

Negative 

electrode 

Separator 

Positive 

electrode 

Thickness  5  [pm] 

50 

25.4 

36.4 

Particle  radius  Rs  [pm] 

1 

- 

1 

Solide  phase  volume 
fraction  solide  £s  [— ] 

0.580 

— 

0.500 

Electrolyte  phase 

volume  fraction  re  [— ] 

0.332 

0.500 

0.330 

Cell  area  Acen  [m2] 

1.0452 

1.0452 

1.0452 

Max  solid  phase 

cone.  csmax  [mol  dm-3] 

16.1 

— 

23.9 

Insertion  rate  cs>e/csmax  at 

SOC  =  0%[-] 

0.216 

— 

0.936 

Insertion  rate  cs>e/csmax 
at  SOC  =  100%  [-] 

0.676 

— 

0.442 

Average  electrolyte 

cone.  Ce,mean  [mol  dm-3] 

1.2 

1.2 

1.2 

Exchange  current 
density  jo  [A  m2] 

36 

— 

26 

Oxydation  transfer 
coefficient  aa  [-] 

0.5 

— 

0.5 

Reduction  transfer 
coefficient  ar  [— ] 

0.5 

— 

0.5 

Solid  phase  Li 

diffusion  coef.  Ds  [m2  s-1] 

2  x  10"16 

— 

3.7  x  10"16 

Electronic  conductivity 
<J s  [S  ITT1] 

100 

— 

10 

Electrolyte  phase 

diffusion  coef.  De  [m2.s_1] 

2.6  x  10"10 

2.6  x  10"10 

2.6  x  10"10 

Electrolyte  ionic 

conductivity  ae  [S  m-1] 

Eq-  (21) 

Eq.  (21) 

Eq.  (21) 

Bruggeman  exponent  p  [— ] 

1.5 

1.5 

1.5 

Li+  transference  number  tp  [— ] 

0.363 

0.363 

0.363 

Table  1  hereafter  details  geometric,  concentration  and  physical 
parameters  required  for  model  simulation.  These  figures  come  from 
the  reference  [4].  Notice  that  the  cell  capacity  is  6  Ah. 

4  A.  Simulating  discharges  at  constant  current 

Fig.  7  presents  cell  voltage  and  insertion  rate  curves  obtained, 
when  simulating  a  complete  discharge  at  6  A  (C  rate),  from  a  fully 
charged  state  (SOC=100%).  Insertion  rates  are  evaluated  on  the 
surface  of  active  material  particles  (r  =  Rs ),  at  electrode/collector 
interfaces  (x  =  0  and  x  =  Lceii).  Curves  in  black  are  associated  with  the 
finite  element  method,  and  curves  in  grey  with  our  model  using 
electrical  analogy.  It  can  be  clearly  observed  that  the  two  models 
nearly  give  the  same  results.  A  really  little  difference  appears  on  cell 
voltages  at  the  end  of  discharge,  mainly  due  to  the  sharp  variation  of 
positive  electrode  OCV  when  insertion  rate  cSie/cs max  is  within  the 
interval  [0.98, 1  ].  For  example,  after  1  h  discharge,  the  finite  element 
method  gives  3.375  V  for  cell  voltage,  0.176  and  0.959  for  insertions 
rates.  Corresponding  figures  obtained  by  electrical  analogy  method 
are  3.376  V,  0.177  and  0.959,  respectively.  Other  example,  a  cell 
voltage  of  3  V  is  obtained  at  t  =  3774  s  with  the  finite  element 
method,  and  at  t  =  3778  s  with  our  analogy  method.  Associated 
insertion  rates  are  nearly  identical.  It  can  be  concluded  that  the  two 
resolution  methods  gives  same  results  for  slow  regimes. 

Fig.  8  presents  the  same  kind  of  comparison  as  previously,  but 
for  a  60  A  (10C  rate)  complete  discharge.  Once  again,  it  can  be 
observed  a  really  fine  agreement  between  the  two  resolution 
methods,  and  a  little  difference  on  cell  voltages  at  the  end  of 
discharge  (a  cell  voltage  of  2.7  V  is  obtained  at  t  =  232  s  with  the 
finite  element  method,  and  at  t  =  234  s  with  our  analogy  method). 
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finite  elem 

ent  method 

electrical  analogy  method 

1  1  1  1  1  1  1  ITH 

Fig.  7.  Simulation  of  a  C  complete  discharge:  cell  voltage  and  insertion  rates  cSie/csmax 
at  electrode/collector  interfaces  (initial  SOC:  100%,  temperature:  300  K). 


Fig.  9.  Simulated  responses  to  a  20C  charge  or  discharge  current  step  (initial  current: 
0  A,  initial  SOC:  43%  (OCV  =  3.6  V),  temperature:  300  I<). 


Last,  Fig.  10  shows  relaxation  curves,  following  a  charge  or 
discharge  regime  at  60  A  (10C  rate)  during  180  s,  the  initial  SOC 
being  0%  or  100%,  respectively.  Once  again,  the  comparison 
between  the  two  models  demonstrates  really  good  agreement. 


4.2.  Simulating  transient  regimes  and  relaxation 

Fig.  9  presents  simulated  voltage  responses  to  a  charge  or 
discharge  current  step  The  initial  current  is  zero,  the  step  amplitude 
is  120  A,  the  initial  SOC  is  43%  (corresponding  to  an  OCV  of  3.6  V), 
and  the  temperature  is  300  K.  It  can  be  observed  a  really  little 
difference  during  the  first  500  ms.  Further  simulations  have  shown 
that  it  was  mainly  due  to  discretization.  Indeed,  results  became 
nearly  identical  if  each  domain  of  the  finite  element  model  is  dis¬ 
cretized  in  10  elements,  as  for  our  model,  instead  of  40  elements. 


4.3.  Advantages  of  the  electrical  analogy  method 

Compared  with  the  finite  element  method,  our  electrical 
analogy  modeling  method  enables  direct  implementation  of  the 
model  in  standard  simulation  software  used  in  electrical  engi¬ 
neering.  Consequently,  it  can  be  easily  employed  for  the  simulation 
electrical  systems  using  lithium-ion  batteries,  as  illustrated  in 
section  5.  One  can  also  take  advantages  of  the  usual  environment  of 
electrical  engineering  simulation  software  to  perform  specific 
analyses,  such  as  parametric  studies,  impedance  spectroscopy, 
investigation  on  charge  or  discharge  strategies. 

Even  in  terms  of  phenomenon  modeling,  our  electrical  analogy 
method  can  provide  facilities.  For  example,  modeling  the  electric 
double  layer  (EDL)  phenomenon  through  a  distributed  capacitive 
coupling  between  solid  and  electrolyte  phases,  in  a  circuit  software 
environment,  is  quite  obvious.  In  our  case,  we  only  need  to  add 
a  capacitance  CDL,eit  between  the  two  circuits  describing  electronic 
and  ionic  charge  transport  in  electrodes  (cf.  Fig.  5),  with: 


Fig.  8.  Simulation  of  a  10C  complete  discharge:  cell  voltage  and  insertion  rates  cs  e/  Fig.  10.  Simulated  cell  voltage  relaxation,  following  a  10C/180  s  charge  or  discharge 
Csmax  at  electrode/collector  interfaces  (initial  SOC:  100%,  temperature:  300  K).  (initial  SOC  before  charge:  0%,  initial  SOC  before  discharge:  100%,  temperature:  300  K). 
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Q)L,elt  =  -y”  CDL>  (35) 

Cdl  being  the  EDL  capacitance  of  the  electrode.  For  illustration, 
Fig.  11  compares  the  simulated  voltage  responses  to  a  20C  charge  or 
discharge  current  step,  as  depicted  in  Fig.  9,  with  the  same  kind  of 
simulation  carried  out  with  description  of  electric  double  layer 
effect.  For  this  purpose,  electrode  EDL  capacitances  are  defined 
versus  the  cell  EDL  capacitance  Cdl, ceil  as  follows: 


Q)L,neg  — 

(i+H 

V  Opos/ 

Q)L,pos  — 

(,+H 

V  °negy 

'  Q)L,cell 
•  Ql, cell 


(36) 


Another  example  is  the  modeling  of  the  solid/electrolyte  inter¬ 
face  (SEI)  layer,  by  an  additional  specific  resistance  Rse\  (in  0  m2) 
which  modifies  the  volume  rate  Li+  current  generation,  as  proposed 
by  Smith  and  Wang  [1,4]: 


ju (x)  =  as  -jo  •  (exp  (aoV^'j  -  exp  (  - 
v(x)  =  ?i(x)  —  ~~*JLi(x) 

This  kind  of  circular  definition  is  not  allowed  by  the  multi¬ 
physics  finite  element  software  we  use,  whereas  it  does  not  present 
any  problem  of  resolution  with  our  electrical  engineering  simula¬ 
tion  software,  which  of  course  includes  a  Newton  iteration 
algorithm. 

5.  Application  on  electrical  power  systems 

In  this  section,  we  give  two  examples  of  use  of  our  lithium-ion 
battery  model  for  simulation  of  electrical  systems.  The  first  aims 
to  show  how  our  model  can  help  battery  conception  (cell  balancing, 
battery  management  system).  It  deals  with  series  association  of 
lithium-ion  cells.  Implementing  such  a  battery  is  really  easy  and 
natural,  when  using  an  electrical  circuit  simulation  software. 
Indeed,  the  single  cell  model  previously  described,  which  appears 
like  a  two-port  component  (in  our  case,  a  current-controlled 
voltage  source,  as  depicted  in  Fig.  6),  can  be  added  to  the  soft¬ 
ware  library,  with  a  symbol  and  a  list  of  parameters.  Then,  it  can  be 
used  in  an  electrical  sketch  like  any  other  component  of  the  library. 


3.64  -i 


3.62- 


<L> 

OXJ 

$  3.6 
o 
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3.58- 


3.56 


-0.5 
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without  EDL  coupling 

with  EDL  coupling 


discharge 


0  0.5 

time  (ms) 


Fig.  11.  Simulated  responses  to  a  20C  charge  or  discharge  current  step  (initial  current: 
0  A,  initial  SOC:  43%  (OCV  =  3.6  V),  cell  EDL  capacitance  :  10  F,  temperature:  300  K). 


In  particular,  it  can  be  associated  in  series  with  other  single  cells, 
with  possible  different  parameter  values.  This  is  what  is  simulated 
hereafter. 

Fig.  12  presents  cell  voltages,  battery  voltage  and  insertion  rate 
curves  obtained,  when  simulating  a  complete  charge  of  a  lithium- 
ion  battery  at  60  A  (10C  rate).  The  battery  consists  of  five 
elements  in  series,  with  identical  initial  SOC.  Solid  phase  diffusion 
coefficients  of  cell  1  is  assumed  to  be  twice  lower  than  those  of 
other  cells,  because  of  a  different  temperature  for  example.  Inser¬ 
tion  rates  are  evaluated  on  the  surface  of  active  material  particles 
(r  =  Rs)t  at  the  electrode/separator  interfaces  (x  =  <5neg  and 
x  =  Lcen  -  <5p0S),  which  are  the  most  critical  places  for  insertion  rate, 
as  far  as  solid  phase  conductivities  (100  S  m-1  at  the  negative 
electrode,  and  10  S  m-1  at  the  positive  electrode)  are  much  greater 
than  electrolyte  ionic  conductivity  (about  0.6  S  m-1).  As  lithium 
diffusion  in  solid  phase  of  cell  1  is  slower  than  in  other  cells,  lithium 
concentration  near  particle  surface  of  cell  1  is  greater,  so  that  the 
voltage  of  this  cell  increases  faster  than  the  others.  As  a  result,  cell  1 
voltage  quickly  reaches  its  rated  value  (about  4  V),  long  before  the 
battery  rated  voltage  (about  20  V)  is  obtained.  Moreover,  in  the 
simulation  carried  out  here,  a  plating  phenomenon  (metal  lithium 
deposition  at  negative  electrode)  can  be  for  seen.  Through  this 
quite  simple  example,  it  is  clear  that  cell  balancing  is  necessary. 

The  second  application,  which  consists  of  instantaneous 
modeling  of  a  power  electronic  converter  involving  a  lithium-ion 
battery,  aims  to  demonstrate  that  our  model  is  suitable  for  elec¬ 
trical  system  simulation.  In  the  case  presented  here,  as  depicted  in 
Fig.  13,  the  converter  is  a  current  regulated  two-quadrant  chopper, 
linking  a  lithium-ion  battery  constituted  with  14  elements  in  series 
(all  elements  are  assumed  identical),  and  a  DC  link  (voltage:  240  V). 
The  input  inductance  value  is  1  mH,  and  switching  frequency  is 
2  kHz.  Battery  current  regulation  is  achieved  by  means  of  a  PI 
corrector,  in  association  with  a  PWM  generator.  Fig.  14  hereafter 
shows  battery  current  and  voltage  waveforms,  in  response  of 
changes  in  current  reference  (0  A  over  [0  s,  5  s],  -60  A  over  [5  s, 
15  s],  0  A  over  [15  s,  20  s],  60  A  over  [20  s,  30  s],  and  0  A  over  [30  s, 
35  s]). 


Fig.  12.  Simulation  of  a  10C  complete  charge  of  a  5-cell  battery  (initial  SOC:  100%,  solid 
phase  diffusion  coefficients  of  cell  1:  Dsneg  m  (2  x  10  16)/2  m2  s_1, 
Ds>Pos  =  (3.7  x  1016)/2  m2  s"1. 
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Fig.  13.  Sketch  of  a  power  electronic  structure  associating  a  14-cell  lithium-ion  battery 
and  a  two-quadrant  chopper. 


Fig.  14.  Simulation  of  an  association  14-cell  lithium-ion  battery/power  electronic 
converter  (initial  SOC:  46%,  temperature:  300  K,  DC  link  voltage:  240  V,  switching 
frequency:  2  kHz). 


Further  work  has  to  be  led  to  complete  the  model  with  an  actual 
experimental  characterization  procedure,  in  order  to  identify 
model  parameters.  The  model  has  also  to  be  studied  versus 
temperature,  and  coupled  with  a  thermal  model  (based  on  elec¬ 
trical  analogy)  of  the  battery,  in  order  to  describe  thermal  runaway 
or  electrothermal  coupling  phenomenon  that  may  arise  under 
some  operating  conditions,  especially  at  low  temperature  (less  than 
0°C). 
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6.  Conclusion 

In  this  article,  we  detail  a  lithium-ion  battery  electrical  model 
that  includes  transient  description,  double-layer  phenomenon,  and 
ID  mathematical  representation  of  charge  and  mass  transport 
phenomena  in  electrodes  and  electrolyte.  For  this  purpose,  elec¬ 
trical  analogies  of  migration  and  diffusion  are  used,  so  that  the 
model  is  directly  implemented  in  simulation  software  dedicated  to 
electrical  circuit  and  electrical  systems  design.  A  comparison  with 
results  obtained  by  means  of  a  finite  element  resolution  method 
have  shown  that  the  use  of  electrical  analogy  for  charge  and  mass 
transport  representation  gives  nearly  equivalent  results.  At  last, 
a  series  operation  and  an  association  with  a  power  electronic 
converter  are  simulated  in  order  to  show  advantages  of  our 
modeling  method  for  system  oriented  simulation.  Of  course,  one  of 
the  main  limit  of  this  method  is  multi-dimension  modeling:  elec¬ 
trical  analogy  is  suitable  for  ID  modeling,  but  not  really  for  2D  or 
3D  modeling. 


Nomenclature 


Latin  letters 

a:  active  surface  per  volume  unit  [m_1] 

A:  active  surface  [m2] 
c:  concentration  [mol  m-3] 

Cdl. '  electric  double  layer  capacitor  [F] 

D:  diffusion  coefficient  of  species  [m2  s-1] 

F:  Faraday’s  constant,  96,485  [C.mol-1] 

GLl:  volume  rate  of  Li+  generation  [mol  m-3  s-1] 
jo:  exchange  current  density  [A  m-2] 
jLl:  volume  rate  of  Li+  current  generation  [A  m-3] 
J :  current  density  [A  m-2] 

I:  current  [A] 

N:  charge  or  species  molar  flow  [mol  s-1] 
p:  Bruggeman  exponent  [— ] 
r:  radial  coordinate  [m] 

R:  radius  [m] 
t:  time  [s] 

tp :  Li+  transference  number  [— ] 

T:  temperature  [K] 
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U:  electrode  equilibrium  potential  [V] 
uf.  thermal  voltage  [V] 

V;  voltage  [V] 

x :  Cartesian  coordinate  [m] 

Greek  letters 

a:  transfer  coefficient  [— ] 

5:  layer  thickness  [m] 
e:  volume  fraction  [— ] 

7]:  surface  overvoltage  [V] 
a:  electrical  conductivity  [S  m-1] 

(p:  charge  or  species  molar  flow  density  [mol  m-2  s-1] 
(f>:  electrical  potential  [V] 


Superscripts  and  subscripts 
cell:  battery  cell 
e:  electrolyte  phase 
eff:  effective  value 
init:  initial  value 
mean:  average  value 
neg:  negative  electrode 
o:  oxydation 
pos:  positive  electrode 
r:  reduction 
s:  solid  phase 

s,e:  solid/electrolyte  interface 
sep:  separator 


